Continuous variable versus EIT-based quantum memories 
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We discuss a general model of a quantum memory for a single light mode in a collective mode 
of atomic oscillators. The model includes interaction Hamiltonians that are of second order in the 
canonical position and momentum operators of the light- and atomic oscillator modes. We also 
consider the possibility of measurement and feedback. We identify an interaction Hamiltonian that 
leads to an ideal mapping by pure unitary evolution and compare several schemes which realize 
this mapping using a common continuous-variable description. In particular we discuss schemes 
based on the off-resonant Faraday effect supplemented by measurement and feedback and proper 
preparation of the atoms in a squeezed state and schemes based on off-resonant Raman coupling as 
well as electromagnetically induced transparency (EIT). 



I. INTRODUCTION 



One of the key elements of quantum communication 
schemes (H HJ and network quantum computers [1, Q 
is a high-fidelity, reversible quantum memory. In recent 
years substantial progress has been made in the design 
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perimental realization 

HI, Hl| of such memories based on ensemble of atoms or 
other quantum radiators as storage unit and photons as 
information carrier. In particular two distinct approaches 
to a quantum light-matter interface have been put for- 
ward: one making use of the off-resonant Faraday effect 
to transfer a qubit encoded in the polarization state of 
a light field to a macroscopic atomic spin of the atomic 
ensemble [1, 1, S B i, [M Gl H, HI , and another one us- 
ing Raman interactions and electromagnetically induced 
transparency, [EIT] lMJI El, El El El El El, H, 
[H, HI [H, |H |13, mmjjr Currently substantial work 
is devoted to the optimization of these schemes which 
is crucial for potential large-scale implementations. An 
obstacle in this effort is the rather different theoretical 
frameworks used to describe these approaches. We here 
put forward a common description of the Faraday and 
Raman schemes, which can be used to compare their ad- 
vantages and drawbacks. 

We begin with a review of general properties of a 
single-mode quantum memory in terms of continuous 
light and matter variables recovering ideal Hamiltonians 
for a purely unitary realization of the map. We then dis- 
cuss the possibility to realize an ideal one-way map from 
light to atoms using Hamiltonians that are not equivalent 
to the ideal one, such as the Faraday interaction Hamilto- 
nian, by means of measurement and feedback techniques 
as well as proper state preparation of the atomic ensemble 
[1, HI] . It is shown that the rather demanding conditions 
on measurement and state preparation can be substan- 
tially reduced in double-pass configurations [II] and can 
be totally eliminated in a triple-pass scheme. 

We then discuss physical implementations of the differ- 



ent mapping approaches starting from the J = — 1/2 <-> 
,/ = +1/2 scheme of Ref. (2(| realizing a Faraday cou- 
pling. We then show that an off-resonant Raman cou- 
pling inaJ = n<->J = ffl configuration allows to 
implement both the Faraday coupling Hamiltonian and 
the ideal mapping Hamiltonian. Finally we discuss the 
quantum memory based on EIT [151 ] and show that this 
scheme corresponds to an ideal mapping Hamiltonian 
when a proper spin polarized (but not squeezed) initial 
state of the atomic ensemble is considered. 

The paper is organized as follows. In Sec. [Til we con- 
sider a simple model for quantum memories in terms of 
atomic and light field quadrature variables. In Sec. Hill 
we discuss physical systems that can be used to realize 
Faraday-type, off-resonant Raman-type, and EIT inter- 
action Hamiltonians. Sec. IIVI summarizes our results. 



II. REALIZATIONS OF THE IDEAL MAP 
USING UNITARY EVOLUTION, 
MEASUREMENT, AND FEEDBACK 

We consider an abstract model of a reversible memory 
for the quantum state of a light mode in an ensemble of 
atoms. The light mode (system L) is described in terms 
of the canonical quadrature variables Xl and P^. We 
assume that the quantum memory (system A) to which 
we intend to transfer the quantum state can be described 
by a similar set of continuous variables Xa and Pa with 
[X^jPa] = i, (% = 1). The latter is the case, e.g., for a 
large ensemble of initially polarized spins, if the excita- 
tion probability of each individual spin is small. The time 
evolution of the two quantum systems can be described 
in the Heisenberg picture by a map that connects the dy- 
namical variables of the systems at some initial time to 
those at a final time t: (X a (0),Pa(0),X l (0),P l (0)) h-> 
(XA(t), PA(t), Xi,(t), Pl(<)). For an ideal quantum mem- 
ory, we require the mapping to be linear in the quadrature 
variables and complete in the sense that the variables of 
one subsystem are mapped only to those of the other. 
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That is, employing the vector notation 

y = (X A ,P A ,X L ,P L ) T , (1) 
the map has the compact form 



y out = My in with M = 



Mi 
M 2 



(2) 



with Mi being 2x2 symplectic, real matrices. The 
matrices need to be symplectic in order to conserve com- 
mutation relations. 

To implement this map we consider a Hamiltonian evo- 
lution that may be supplemented by measurement and 
feedback processes. To ensure the linearity of the map, 
the Hamiltonian should be of at most second order in 
the quadrature variables (or in the corresponding anni- 
hilation and creation operators). Specifically, we consider 
pure harmonic oscillators with quadratic interaction be- 
tween them: 



H = H A + H L +H mt 
Ha 



(3) 



~"X 2 A + P\) = u A [a\a A + ^), (4) 



Hl = -f[Xl+Pi)=Lj L \a\a L + ^J, (5) 
Hint = pX A X L + qX A P L + rP A X L + sP A P L , (6) 



1 



where lu a and ljl are the oscillator frequencies and the 
real parameters p, q, r, and s characterize the interaction. 
These parameters may have explicit time dependence. 



A. Purely unitary evolution 

Let us first consider the question under what condi- 
tions an ideal quantum memory map can be realized by 
pure unitary evolution [30l [3l| . For this we assume that 
the free Hamiltonian Hq = H A + Hl commutes with the 
interaction Hamiltonian, 



[H (t),H mt (t)]=0. 



(7) 



With this restriction, the free Hamiltonian can be elim- 
inated from the equations of motion. Indeed, let us 
express operators of the Heisenberg picture in 
the frame rotating according to the free Hamiltonian: 
A(t) = U (t)AH(t)Ul(t), where Uo(t) denotes the unitary 
operator of the interaction-free time evolution. Thus, the 
equation of motion for operators in the rotating frame 
reads 



-A(t)=i[H int (t),A(t)]. 



(8) 



If we further assume that the interaction Hamiltonian 
commutes at different times, [H, nt (t) , Hi nt (t')] = 0, then 
the time evolution operator corresponding to the interac- 
tion can be written in an exponential form with no time 



ordering necessary. This assumption together with ([7]) 
necessarily implies exact resonance between the atomic 
and light systems (lo a — ljl — lj) and that the interac- 
tion Hamiltonian is of the form H mt (t) = a(t)Hi with 

H x = S mZ(X A X L +P A P L ) + cosZ(P A X L -X A P L ), (9) 

that is, i?int can have explicit time dependence only 
through a(t). Then Eq. ([5]) can formally be solved using 
the Baker-Campbell-Hausdorff formula 



a(r) 



(10) 



where A n = (—i) n [[[A(0), Hi], Hi], . . . ] is proportional to 
the n-fold commutator of A(0) with Hi. It is easy to see 
that the commutators of the quadrature variables y with 
the interaction Hamiltonian © are linear in the same set 
of quadratures, namely, we have 



y,Hi 

where 
C = 



= iCy and y„ = C y n ^ = C n y, (11) 




R 



R 1 




and R 



cos £ — sin £ 
sin £ cos £ 



(12) 

Using the series (fTO]) and recognizing that C 2 = — 1, we 
find the linear relation between the quadratures at time 
t and the initial time 



m 

with 



cos$(£) 1 + sin$(t) 



-R- 1 
R 



dr a(r). 



y(o), (13) 



(14) 



We see that perfect quantum memory mapping is 
achieved after an interaction time T such that the area 
of coupling is $(T) = (2n + l)vr/2, (with n € Z). In this 
case one has 



y(T) = ± 



-R- 1 
R 



y(o). 



(15) 



Regarding the possibility of physical realizations of the 
interaction Hamiltonian ^ , two different values of £ are 
of particular importance. For £ = 0, we have 



H int (t)=a(t) (P a Xl-X a P l 



(16) 



If the envelope a(t) is chosen such that $(T) = tt/2 then 
we arrive at an ideal quantum memory map 



X A (T)=X L , P A {T)=P Ll 
X L (T) = -X A , P L {T) = -P A . 



(17) 
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We can also consider the storing process in which the 
X and P quadratures of the memory system are in- 
terchanged with respect to the previous transformation. 
This corresponds to £ = 7r/2 and the interaction Hamil- 
tonian 

flint (*) = a(t) (x A X L + P A P L ^ (18) 
leads to the map 

X A (T)=P L , P A (T) = -X L , 
X L (T) = P A , P L (T) = -X A , 

which is again the map of an ideal quantum memory. 

B. Single-pass scheme with feedback and initial 
spin squeezing 

A quantum memory for light, which is not based en- 
tirely on unitary evolution but is rather an approximate 
simulation of the Hamiltonian (|18l) . was proposed and 
experimentally demonstrated in [201 ] . The light-matter 
interaction there is due to the Faraday effect and is de- 
scribed by the Hamiltonian 

H = const., ffi = P A P L . (20) 

See Sec. MI Al for a possible derivation of (|2H)) . In this 
case the matrix C cannot be represented in the form (TT2|) 
and the unitary evolution will simply shift the position 
operators by an amount proportional to the momentum 
of the other system as well as to interaction time t, while 
the momenta are constants of motion: 

X A {t)=X A +tP L , P A {t)=P A , 
X L (t)=X L +tP A , P L {t) = P L . ( " 1 

One recognizes that only the momentum quadrature of 
the light mode is transferred to the atomic ensemble. To 
also map the position quadrature to the ensemble the uni- 
tary evolution was complemented in Ref. 20] by a homo- 
dyne measurement of the outgoing light quadrature Xl 
and the measurement result x was fed back by applying 
a momentum displacement of —x/t on system A. As a 
consequence of the measurement, one can formally write 
the c-number x in place of X^lt) and rearrange (|2ip to 
conclude that after measurement and feed-back: 

XT m — X A + tP L , Pr m = -\*l- ( 22 ) 

If the atomic ensemble were initially prepared in a posi- 
tion eigenstate, i.e., in an infinitely squeezed state, the 
operator X A in (f2"2")l could be replaced by a c-number. 
The resulting map would in this case ideally transfer the 
complete state of light to the atomic ensemble. 

To verify these statements in a more rigorous way, we 
calculate the state of the quantum memory after the stor- 
age and the storage fidelity in terms of Wigner functions. 



The atomic and light systems are initially disentangled, 
so the two-particle Wigner function is of the product form 

W (x A)PA ;x L ,p L ) = W A {x A ,p A )W in (x L ,p L ), (23) 

After the transformation (f2Tj) . the new state is given by 

W(x A ,p A ;x L ,p L ) = W (x A - tp L ,p A ;x L - tp A ,p L ). 

(24) 

Then the quadrature X' L of the outgoing light is mea- 
sured. For an ideal measurement, the projection cor- 
responding to the outcome x is 1a ® T\ x , where II^ = 
\x)ll(x\ and its Wigner function is H x (x' L ,p' L ) — S(x' L — 
x). The (unnormalized) conditional atomic state can 
be obtained by a von Neumann projection. The feed- 
back is described by a shift in the atomic momentum 
p'a = Pa + x ft ano - thus the state of the memory condi- 
tioned on the measurement result x and after feedback 
reads: 

WT^/a) = / dx' L dp' L W A (x' A - tp' L ,p" A + x/t) 

x W in (x' L - tp" A - x,p' L )U x (x' L ,p' L ). (25) 

The norm of (|25|) gives the probability distribution 
P(x) of the measurement outcome x. In general, P(x) 
has explicit dependence on the unknown input state 
W 1ti {xl:Pl)- This means that the measurement yields 
information about the input state of the light which is, 
therefore, distorted. However, in the special case of an 
infinitely squeezed initial atomic state W A (x A ,p A ) — 
8{x A — Xo), we have a uniform probability distribution 
P(x) = 1/t and no information about the initial light 
state is obtained by the measurement. In this ideal case, 
the final state is 

W mcm (x,p) = W in {~tp, (x - x )/t) (26) 

that corresponds to an ideal quantum memory mapping 
regardless of the outcome x. 

In reality, however, it is impossible to prepare the 
atomic ensemble in a position eigenstate, i.e., in an in- 
finitely squeezed state. Furthermore, the measurement 
will be imperfect in general. In the followings, we discuss 
the effects of finite initial squeezing and finite detection 
efficiency. For this we replace the Wigner function H x 
used above with the corresponding expression for a finite 
detection efficiency [32[ 

fl x , av = (2^l)- 1 ' 2 J dye~^\y) LL {yl 
^,a v (x' L ,p' L ) = (2^)- 1 /2 exp |_ (<^) 2 |, (27) 

Here a n characterizes the resolution of the position mea- 
surement, with <r„ = corresponding to a noiseless, per- 
fect measurement, while a typical experimental value is 
in the order of a'i = 2.5 x 10~ 3 . Averaging over the 
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measurement outcome x gives the memory state 

W mcm (x,p) = J dx'dp'dx"W A (x',p') 

x W in {-tp + a v x" , (x - x')/t)H 0A (x" ,0). (28) 

Let us consider as the initial atomic state a Gaussian 
spin squeezed state 

-i _ (^-^o) 2 _ (p-po) 2 

W A (x,p) = - e 2 ^ ^ . (29) 

In the experiment of [20| , it was actually a coherent spin 
state with o 2 x A — a 2 P A = i . The imperfect measure- 
ment (|2~T1) then results the state 

i /• _,' ! V 

FF mcm (x,p) = / dx'dp'e ^ 

2na ri ax,A J 

xW in {-tp + x',{x-x )/t + p') (30) 

for the atomic quantum memory. Compared to the ideal 
memory state (|26p . we have higher uncertainties both in 
the X and P quadratures due to the noisy measurement 
and the imperfect initial state preparation, respectively. 
The average fidelity of the storage process can be cal- 
culated from the overlap of the ideal (f2"6")l and real ([ST)]) 
output states, 

F = 2tt J dxdpW in (-tp,(x - x )/t)W mcm (x,p). (31) 

For a Gaussian state of the input light field 

, _(x-x in ) 2 _ (p-p in ) 2 

W in (x,p) = - e 2CT ^ , (32) 

the fidelity reads 

F = [(* x , A /t 2 + 2o* L )(o* + 2a xx )}-^, (33) 

which is about 82% for a coherent spin state and coherent 
light input {<J X a = a XL = °p l ~ V^, o - ?? = 0, and 
t = 1), as in Ref. [20(. However, direct calculation shows 
that the fidelity of storage quickly decreases for highly 
nonclassical states like Schrodinger cat states (see Fig.Q]). 

We conclude that the single-pass feedback technique 
provides the perfect quantum memory mapping if and 
only if the atomic ensemble is initially prepared in an 
infinitely squeezed spin state W A {x A ,p A ) = 5{xa — %o) 
and if the measurement is noiseless. 

C. Double- pass schemes 

The necessity of preparing the atomic ensemble in a 
highly squeezed state or performing a measurement with 
feedback can be avoided in a double-pass scheme with 
two successive, different unitary evolution (To| . In this 




FIG. 1: Fidelity of storing the odd Schrodinger cat states 
\a,— ) oc \a) — \—a) (with real a) in the one-pass scheme 
of [20]. The atomic ensemble is initially prepared in a spin 
squeezed state: the smaller the position variance ox, a, the 
higher the squeezing. Thick line corresponds to the coherent 
spin state with no squeezing and it reaches the classical limit 
of fidelity 0.5 (thick dashed contour line) at about \a\ 2 = 2.0. 
Since the measurement alters the input state, the scheme can- 
not efficiently store superpositions of states with very different 
momenta. 

scheme, an interaction Hamiltonian identical to (120[) is 
applied first. After a time t, the interaction is suddenly 
changed to H 2 , so we have 

Hi = P a Pl, H 2 = X A X L , (34) 
X A {t + t') = X A +tP L , 

PA{t + t , ) = {l-tt l )P A -t'X Ll 

(35) 

X L (t + t')=X L +tP A , 

P L {t + t r ) = (l-tt r )P L -t'XA. 

If the interaction times are adjusted such that tt' = 1, 
we directly obtain a mapping like (|22p . however, without 
measurement and feedback, 

X' A = X A + tP L , P A = -jX L , 

X' L = X L + tP A , P' L = -jX A . (36) 

Perfect mapping can thus be achieved if the atomic 
ensemble is initially prepared in an infinitely squeezed 
state. In such a scheme there is no need for measure- 
ment and feedback, and the fidelity of the memory is 
F = (cr X A /t 2 + l)- 1 / 2 for coherent input light. 

Alternatively, as will be shown in the following, the 
necessity of an initial atomic squeezing can be avoided 
by using measurement and feedback. Indeed, a measure- 
ment of the quadrature P' L can project the initial atomic 
state to a squeezed state — thus appropriately performing 
the atomic state preparation after the interaction [2l| . If 
the measurement of P' L gives a value p, we can replace 
X A by —tp. Applying a position displacement of tp of 
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the atomic position, X^ cm = X' A + tp, we obtain an 
ideal quantum memory map 



XT m = tP L , 



P>mcm 
A 



-Xr 



t 



(37) 



Note that an imperfect measurement similar to (|27[) also 
introduces noise in the position quadrature, 



W mcm (x,p) 



1 



dx'W in {-tp,x/t + x')e 



(38) 

However, light measurement can be performed with far 
much higher accuracy than spin squeezing. For coherent 
input light, the fidelity of storage reads F = (cr 2 + 1) -1 / 2 . 

Thus in the double-pass scheme one can get rid of ei- 
ther the measurement and feedback or the preparation 
of the atomic ensemble in a squeezed state. 



D. Triple- pass scheme 

Finally, we mention that the ideal unitary evolution 
given in Sec. Ill Al can be equivalently achieved in a three- 
pass scheme without measurement and with no initial 
atomic squeezing. As it was pointed out in (3lj . a beam 
splitter like interaction Hamiltonian can be simulated by 
successively applying the above two kinds of Hamiltoni- 
ans three times. In fact if we reapply H\ for a third time, 
we obtain 

H 1= P A P L , H 2 = X A X Ll H 3 = P A P Ll (39) 
X A (t + t' + t") = (1 - t't")X A + [t + i"(l - tt')]P L , 

P A (t + t' + t") = (1 - tt')P A - t'X L , 

X L (t + t' + t") = (1 - t't")X L + [t + t"{\ - tt')]P A , 

P L {t + t' + 1") = (1 - tt')P L - t'X A . 

(40) 

Setting t = t' = t" = 1 we arrive at the ideal mapping 
(TT9")) . Let us denote with U the unitary operator describ- 
ing the time evolution due to the consecutive actions of 
Hi, H 2 , and H 3 over periods of time t — 1, t' — 1, and 
t" = 1, respectively, 

U = e -iPAPt e -i^AXt e -iPAP Lj ( 41 ) 

and let us denote with J7ideai the unitary evolution corre- 
sponding to the Hamiltonian of the ideal quantum mem- 
ory (|18[) with a(t) = 1 over a period of T = ir/2, 



(42) 



The two unitary operators have the same effect on the 
quadrature variables and U and C/jdeai are in fact identi- 
cal, 





* / \ |l) N A /2 N A /2 |2> 
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FIG. 2: Quantum memory scheme based on Faraday rotation. 
The atoms with doubly degenerate ground states are polarized 
in the polarization direction of the strong classical light field. 
The inset shows a simplified atomic level structure when the 
quantization axis is the direction of light propagation. The 
atoms are in an equal-weighted coherent superposition of the 
ground levels. The presence of a weak y-polarized light field 
results in a imbalance of the dynamic Stark shifts for the two 
levels, developing a relative phase in the superposition. The 
atomic spins are, thus, coherently rotated in the xy plane as 
a back-action of light on atoms. 



III. PHYSICAL SYSTEMS REALIZING 
ATOMIC QUANTUM MEMORY 

In this section, we analyze three kinds of configura- 
tions that can serve as collective atomic quantum mem- 
ory for light. The first of them is experimentally car- 
ried out by Julsgaard et al. [2(J. Then we discuss a 
scheme based on off-resonant Raman scattering that re- 
alizes the ideal interaction Hamiltonian (fl8|) directly. Fi- 
nally, we show that electromagnetically induced trans- 
parency (EIT) gives rise to an effective interaction Hamil- 
tonian of type |H 



,-iP A PL p -iX A X L p -%P A Ph - p -i%{X A X L +P A P L ) 



(43) 



A. Quantum memory based on Faraday rotation 

Let us analyze first a scheme in which the light-matter 
interaction originates in the paramagnetic Faraday ef- 
fect [33l ]: given an ensemble of atoms with macroscopic 
magnetic moment and shined by a linearly polarized light 
beam propagating in the direction of the magnetic mo- 
ment, the plane of light polarization is rotated. 

To describe this atom-light interaction, consider an 
atomic level structure depicted in the inset of Fig. [2] 
where the two ground levels are off-resonantly coupled 
to the upper ones by the right and left circularly polar- 
ized electromagnetic field modes of the same frequencies 
wjr = lol = u). In the experiment of Julsgaard et al. 
[20], levels |1) and |2) correspond to the Zeeman sub- 
levels M F = ±4 of the ground state 6 2 S*i/ 2 (F = 4) of 
cesium, while the light pulses are detuned to the blue 
by 700 MHz from the 6 2 S 1/2 (F = 4) -> 6 2 P 3/2 (F = 5) 
transition (A = 852 nm). 
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If the detuning A from the atomic transitions is large 
enough (g(a'a) <C A), we can adiabatically eliminate 
the off-resonant excited levels so that the system reduces 
to an effective two-level atom. The dynamics is then 
governed by the following effective Hamiltonian, 

H = u)(a' R a R + a\a L + I) 

N A 

+£0^(11)^(11 + 12)«<2|), (44a) 

i=l 
N A 2 

= ^|;[4 a fll 1 )»< 1 l + 4^l 2 )n(2|]. (44b) 

where Na is the number of atoms and the electric dipole 
coupling constant for the single-photon transitions is 

g = yZj(l|de B |4)| = ^L\(2\de L \3)\. (45) 

Note that (|44b[) is nothing more than the dynamic Stark 
shift (light shift) of the energy of the lower levels caused 
by virtual transitions to the off-resonant upper levels. 

Let us introduce the quantum mechanical Stokes pa- 
rameters to describe the polarization state of light, 

S x = (a R a L + a' L a R )/2 = (a x a x - a y a v )/2, 

S y = (a R a L - a\a R )/{2i) = {a) x a y + a\a x )/2, 

S z = (a R a R - a\a L )/2 = {a ] x a y - a\a x )/(2i), 

N s = (a R a R + a\a L )/2 = {a) x a x + a\a y ) /2. (46) 

S x is the photon number (intensity) difference of the x 
and y linearly polarized light components, S y is that of 
the diagonally polarized ones, S z corresponds to the dif- 
ference in the right and left circularly polarized com- 
ponents, while Ng is half of the total photon number 
in the two modes. It is easy to see that they satisfy 
the standard angular momentum commutation relations 
[Sa^Sfj] — ie a /3jS 7 and [JVs,S] = 0. The annihilation 
operators a x — (a R + a^)/\/2 and a y — (a R — a,L)/(iy/2) 
correspond to the linearly polarized components of the 
light beam. Note that S is not a vector operator, i.e., it 
does not transform as a vector under rotations. 

For the atomic magnetic moment, we can use 
Schwinger's representation and introduce the collective 
atomic quasi-spin variables, 

u x ee £(|2>«<l| + |l)«<2|)/2, 

i 

&v = EO 2 )^ 1 !-!!)^!)^, 

i 

a z EE ^(|2)«<2|-|l)«(l|)/2, 

i 

N a ee 2(|2)«<2| + |l)«(l|)/2, (47) 



where 2N„ is equal to the number of atoms in the sub- 
space spanned by |1) and |2), i.e., N a counts the atoms 
contributing to the phenomenon. Note that unless we 
have a real spin system, the quasi-spin vector is not a 
vector operator either because [a a , Fp\ ^ ie a p 1 a 1 with 
F being the total angular momentum. By definition, the 
components satisfy the commutation relations 

[o'a, &p] = ieapy&y. (48) 

We are now ready to express the effective interaction 
Hamiltonian (|44b[) in terms of the Stokes and quasi-spin 
vectors, 

H = u(2N s + l) + 2E Na, (49a) 
H int = \ (NsN. - S z a z y (49b) 

The operators Ns, N a , and NsN a all commute with both 
S and (j (and thus, also with S z a z ) and they have nothing 
to do with the dynamics. Therefore, we can safely omit 
them and write 

2g 2 - 

H = const., H int = — -£-S z a z . (50) 

The remaining term explains paramagnetic Faraday ro- 
tation [33j . If the z component of the collective atomic 
quasi-spin has a macroscopic expectation value, it results 
in an interaction Hamiltonian proportional to S z which 
in turn introduces rotation of the Stokes vector along 
the z axis — thus turning the linear polarization in the xy 
plane. Reversely, as a back-action of light on atoms, if 
the z component of the Stokes vector has a macroscopic 
expectation value, it will rotate coherently each atomic 
spin along the z axis. 

Now we show — following Ref. [2(| — that Eq. (f5T7|) gives 
rise to a Hamiltonian of the form (|2H)l . Consider an 
atomic ensemble in which all atoms are initially prepared 
in the superposition state (| 1) + \2))/y/2. In this coherent 
spin state (CSS), the x component of the quasi-spin vec- 
tor has a macroscopic expectation value {<t x )o = Na/2. 
This state is actually a quasi-spin--^ 4 state, i.e., an eigen- 
state of or 2 with eigenvalue ■^■(■^ 1 + 1). Since the Hamil- 
tonian (j49|) commutes with <r 2 , the atomic state will al- 
ways remain a quasi-spin-^ 4 state. Moreover, as long 
as the interaction introduces small perturbation to <r, 
the y and z components of <r stay small with respect to 
the x component. Therefore, we can use the Holstein- 
Primakoff approximation [34| and introduce the position 
and momentum-like quadrature operators 

X A ee 6r y J yj (a x ) , Pa = o z j \/ (a x ) . (51) 

As long as the number of atomic excitations ha — 
h(X^ + P| — 1) is small compared to the number of 
atoms, the deviation of the quasi-spin vector from the 
CSS stays in the tangent plane of the spin sphere. Then 
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Eq. (|48)l ensures that we have the approximately correct 
commutation relation 



[X a ,Pa] 



2i 
Na~ 



n A 



(52) 



In the experiment of [20(, the quantum information is 
represented in the y-polarized weak signal beam propa- 
gating in the z direction and it is mixed with the coprop- 
agating x-polarized strong classical control field having 
coherent amplitude a (Fig. Writing the c- number 
a in place of a x , we find that the x component of the 
Stokes vector (|4"6")l has a macroscopic expectation value, 
(S x ) = |a| 2 /2 in first order. Therefore, we can introduce 
the light quadrature variables as 



X L = S y /J(S x ) = (e-^ay + e^&D/Vz.. 



Pt. 



S z / /S)o = (e-*% - e**&l)/iy/2, (53) 



where ip = arg(a) is the complex phase of a. 

In summary we express the Stokes and quasi-spin vec- 
tors in terms of the quadrature operators [cf. Eqs 
and (jIT)) ]. 



S x = \a\ 2 /2-h Ll 



N A /2 - n A , 



S y = 



& y = XaVNa/2, 



N s = H 2 /2, 



a z = PaVNa/2, 
N a = N A /2. 



(54) 



Note that these operators exactly satisfy the angu- 
lar momentum-like commutation relations [-/V ct ,<t] = 0, 

[<ra,vp] = jea/3 7 ^7, [Ns, S] = 0, and [S Q ,S/3] = ie a /3 7 S^, 
despite the fact that we used an oscillator approxima- 
tion for the atomic system and a classical model for the 
x-polarized light. 

Putting ([51]) and ([53]) into ([50)) we arrive at the in- 
teraction Hamiltonian in terms of the atomic and light 
quadrature operators, 



Hint — — : 



A 



-PlPa. 



(55) 



We recover the non-ideal interaction Hamiltonian (|2H)l . 
As a consequence, the atomic position quadrature Xa is 
displaced by an amount proportional to Pl [cf. Eq. (|2T|) ]. 
thus rotating the quasi-spin vector towards the y axis. 



B. Off-resonant Raman scheme 

Now we discuss another physical system which directly 
realizes the ideal interaction Hamiltonian (fT8|) . Our 
quantum memory scheme is based on off-resonant Ra- 
man scattering by A- type atoms 

Consider an ensemble of Na atoms depicted in Fig. [3] 
that is shined by a monochromatic light beam consisting 



right polarized 
classical field 




FIG. 3: Quantum memory scheme based on off-resonant Ra- 
man scattering. The inset shows the A-type atomic config- 
uration. The degenerate ground states |1) and |2) are off- 
resonantly coupled to the intermediate levels |3) and |4) via 
right and left circularly polarized monochromatic light beams. 
Light shifts of the lower levels are cancelled by tuning the 
laser fields right between the two upper levels. The classi- 
cal control field is the right polarized one. Only level |2) is 
macroscopically populated, which corresponds to a quasi-spin 
polarization in the direction of light propagation. Absorption 
of a signal photon results a collective atomic excitation in 
level |1). 



of copropagating phase-locked right and left circularly 
polarized components. Levels |1) and |2) can be, for ex- 
ample, the itif = ±1 Zeeman sublevels of the 2 5i/2 (F = 
1) ground state of sodium or rubidium (with nuclear spin 
I = 3/2), while levels |3) and |4) are the mp = Zee- 
man sublevels of the hyperfine levels 2 P\/2{F = 1) and 
2 P\/2 (F = 2). Although the Zeeman sublevels M F = ±2 
of the uppermost level 2 Pi/2 (F — 2) also contributes to 
the light shifts, the principle of the model is not changed. 

If the atomic transitions are far off-resonant, the four- 
level atom is reduced to an effective two-level one. Af- 
ter adiabatic elimination of the upper levels, the dynam- 
ics is governed by an effective interaction Hamiltonian 
that consists of dynamic Stark shifts and two-photon pro- 
cesses between the ground states, 



N A ( 

1=1 



(\9m? 



A 

A 12 



+ 1— + ^J4^|2>»<2| 

/ 93B.93L . 9ar9al Vt - io\ /1 1 . xi „ 
I ^ + — 2y — Ja' L a R \2} u {l\ +H.c. 



, (56) 



where the electric dipole coupling constants for the 
single-photon transitions are 



93R 



2HV 



;i|d-e R |3) e 



i(k R ri+0 H ) 



(57) 



and similarly for g l 3Ll g\ R , and g\ L , with being the 
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position of the ith atom. For symmetry reasons, we 
have \gi R \ = \g\ L \ = g and \g\ R \ = \g\ L \ = g'. Fur- 
thermore, since the frequencies and propagation direc- 
tions of the two polarized light beams coincide, we can 
choose the relative phase 4>R ~ 4>l of the mode functions 
of the light modes so that g\ R g\^* — g 2 is real and pos- 
itive. However, in our example of sodium or rubidium, 
g\ R and g\ L have opposite signs for this choice of phases, 

so 9 1 4r9\l* = -a 12 - 

We can introduce the photonic Stokes vector S PU)) 
and the collective atomic quasi-spin vector a (|47|) . If the 
ground states are F = 1 hyperfine sublevels, the com- 
ponents of the quasi-spin expressed by the total angular 
momentum F are 



i 



a z = F z /2. 



(58) 



We are now ready to simplify the interaction Hamilto- 
nian ()56[> . We realize that the first two terms therein are 
proportional to NsN a — S z a z and they are responsible 
for paramagnetic Faraday rotation. Since NgN a com- 
mutes with both S and a and has nothing to do with the 
dynamics, we will omit it in the following. The last term 
in (I56p is proportional to S x a x + S y a y . All in all, we have 



f£_ 
A 



A' 



<£_ 
A 



S z a z 

_ 9 
A 



-}j (S x a x 



Sy&y 



(59) 



We note that the effective Hamiltonian is essentially of 
the same form even for more gen eral configurations in- 
volving multiple atomic levels [351 ]. The meaning of the 
first term is explained in Sec. IIII Al By tuning the 
laser fields right between the two upper levels so that 
g' 2 1 A' — —g 2 /A, we can cancel the first term. The re- 
maining term then describes the two-photon processes of 
the atoms making transitions between the two ground 
states. However, this term can also be interpreted as 
Faraday rotation. Suppose, for example, that the y com- 
ponent of the Stokes vector has a macroscopic expecta- 
tion value, that is, the atomic ensemble is irradiated by 
a single 45° linear polarized classical beam. Then the 
interaction Hamiltonian (|59p reduces to a term propor- 
tional to <j y corresponding to a coherent rotation of the 
quasi-spin along the y axis. Indeed, if we choose y as the 
quantization axis, such a diagonally polarized light will 
induce virtual atomic transitions from the superposition 
state (|l)+i|2))/v2 to the off-resonant upper levels, thus 
shifting the energy level of this superposition state with 
respect to the orthogonal state (|1) — i|2))/V2- As a 
consequence, all the atomic quasi-spins are rotated along 
the y-axis. Appropriately choosing the polarization and 



detuning of a single classical light field, one can realize 
rotation of the quasi-spin vector along arbitrary axis. 

To use the present configuration as a quantum mem- 
ory, the polarizations of the control and signal fields and 
the initial atomic state should be chosen differently from 
that of Sec. IIII Al Let us represent the quantum infor- 
mation in the weak left circularly polarized light beam 
(signal) and let the right circularly polarized light be the 
strong classical field (control) with coherent amplitude a 
(though not too strong so that it remains off-resonant, 
\a\ <C |A/g|, \A'/g'\). We can write the c-number a 
instead of the right circularly polarized annihilation op- 
erator a R and we find that the z component of the Stokes 
vector is a classical variable (S z )o = M 2 /2- This enables 
us to introduce the light quadrature variables as 



X L = S x /JW/2 = (e-^a L + e^a{)/V2, 

P L = 5„/v/HV2= {e-^a L ~e rv a{)/iV2, (60) 

where tp = arg(a) is the phase of a. Similarly we can 
introduce quadratures for the atomic ensemble if the col- 
lective atomic state stays close to the coherent spin state 
(CSS) in which all atoms are in level \2). The z com- 
ponent of the quasi-spin has a macroscopic expectation 
value, (a z )o = Na/2 in zeroth order, and we can approx- 
imate it with a z — Na/2 — fiA- Then we can write the 
atomic quadrature variables as 



X A = 



VNl/2 



Pa = 



VNa/2 



(61) 



In order to keep the correct angular momentum com- 
mutation relations of the Stokes and quasi-spin vectors, 
we write 



S x =X L ^\a\ 2 /2, 



= XaVnI/2, 



Sy=P L VW/2, 

S z = \a\ 2 /2~n L , 
N s - |a| 2 /2, 



a y = PaVNa/2, 
Oz = Na/2 - fiA, 
N a = Na/2. 



(62) 



Note that the number of atomic excitations ha = + 

P\ — 1) equals to the population of level |1), ua = 
Putting all together while canceling the term 
S z g z in J59|) . we obtain the interaction Hamiltonian 



2g 2 \a\ 



-{X l Xa + PlPa)- 



(63) 



The Hamiltonian (|63|) results in an oscillation of the 
excitations between modes A and L. Indeed, with the 
creation and annihilation operators of systems A and L 
we recognize the beam splitter Hamiltonian 



X L X A + PlPa = a\d A 



a A a L - 



(64) 



Absorption of a left polarized photon makes an atomic 
transition from level \2) to level |1), and reversely, emis- 
sion of such a photon causes a decrease in the population 
fiA of level |1). 
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13) 




A 



n A 




N A - n A 



FIG. 4: A-type atomic configuration in the EIT scheme. Lev- 
els |1) and 1 3) are coupled by the classical field of Rabi fre- 
quency Q(t) which is controlled externally. Levels |2) and |3) 
are coupled by the weak signal field. Initially, only level |2) is 
populated. 



The state of the two systems can be exchanged com- 
pletely by appropriately adjusting the amplitude a(t) of 
the classical right circularly polarized control field, so 
that we have a w/2 pulse, 



/ \a(t)\dt = -. 



(65) 



This leads to the ideal quantum memory mapping (fT9 



C. Resonant EIT scheme 

In quantum memories based on EIT [l5| . an intense 
classical radiation field (control) and a weak quantum 
field to be stored (signal) are adjusted on or near res- 
onance with the transitions of the A-type atoms (see 
Fig. 3]). In this subsection we show how to adiabati- 
cally eliminate the resonant excited level and we derive a 
beam splitter-like effective Hamiltonian between the sig- 
nal light field and the collective coherences of the lower 
atomic levels. In order to simplify the discussion, we will 
restrict ourselves to single-mode radiation fields for the 
control and signal with exact two-photon resonance and 
standing polariton wave. The Hamiltonian in the rotat- 
ing frame then reads 



H 



N A 

E 

i=l 



[-A|3)«(3| + ( Si 6£|3)«<2| + n < |3) < 



(ll+H.c. 



(66) 

where is the bosonic operator of the signal field, A is 
the one-photon detuning for both transitions. Further- 
more, we disregard atomic motion and assume that the 
coupling constants <?i are real and the same for all atoms, 
and so are the Rabi frequencies fli. This allows us to in- 
troduce collective spin operators 

<7 o6 = ^|a>«<6| (a, b =1,2, 3). (67) 

i 

Due to the symmetry of the Hamiltonian (f6"6")l , only the 
totally symmetric Dicke states are coupled to the light 



fields. The totally symmetric state containing n atoms 
on level |1), m atoms on level |2), and I on level |3) is 
defined as 



\ni, 77^2 1 ^3) = \n\ m\ l\ (n + m + 



-1/2 



x |1: ii, . . . , 3: k u ...,ki), (68) 

where the summation is over the N A = n+m+l mutually 
different variables h, . . . , i n , ji, . . . , j m , hi, . . . , ki going 
from 1 to N A . The ket vector at the right of (f6"5)) repre- 
sents the atomic product state in which atoms indexed 
by «i, . . . , i n are in state |1), etc. It is easy to verify that 
the spin-flip operators l|67p keep the symmetry of these 
states. 



a ab \n a ,m b ,l c ) = y/ {n + l)m\(n + l) a , (m - l)f,, l c ), 
Vaa\n a ,m b ,l c ) = n\n a ,m b ,l c ), (69) 

with {a, b, c} = {1, 2, 3}. We can see that for each value 
of n = 0, . . . , N A , the set 

S n = Span {\(k - 0i, (N A - k) 2 , l 3 ) \n - k) L \ 

k = 0, . . . , min(n, N A ) and I = 0, . . . , k) (70) 



is an invariant subspace of the Hamiltonian (|66|) so that 
([6^| is block diagonal. (The subscript L means that the 
corresponding Fock state refers to the signal light mode.) 

Now we make use of the fact that only the atomic 
level 1 2) is macroscopically populated, so 022 can be sub- 
stituted by the c- number N A . We observe that only 
the coherences CT21 and 023 are important, as we can 
express the other spin operators as &n — A , 
CT33 = 0320-23/NA, and CT13 = CTi2(T23/^A- Compared 
to the single atomic oscillator in the off-resonant Ra- 
man process investigated in Sec. IIII B[ now we have 
two atomic oscillator modes with annihilation operators 
<J2i/\/N A and <J23/VN A , respectively. Within the sub- 
space S = © n< jVA &n of photonic and totally symmetric 
atomic states with arbitrary but small number of exci- 
tations 77, the Hamiltonian can be written in the matrix 
form fal 



H 




(71) 



Here we have introduced the vector notation = 
&i2/\/N A , (732/V-^A) for the creation operators of 
the photonic mode and the two atomic oscillator modes. 
Note that in the limit of small atomic excitation, compo- 
nents of z and z^ approximately satisfy bosonic commu- 
tation relations. The matrix H can be brought to a block 
diagonal form. The corresponding transformation matrix 
R defines the quantum field variables (annihilation oper- 
ators) of the so-called dark- and bright-state polaritons 

'cos 6 -sin6» 0^ 
Rz, R = ( sin6> cos<9 | , (72) 
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where $, and S stand for the dark-polariton, the 
bright-polariton, and excited-state modes, respectively. 
The mixing angle 9 is defined as 



tan 9 



(73) 



Note that in the limit of small atomic excitations, the 
polariton operators also satisfy the bosonic commutation 
relations, 

[*,^+] « « |^,Ht] » 1, 

[*,*t]„ [$,Ht]« [$,Ht]«0, 
[*,$] = [*,H] = [*,§]=0. (74) 

Actually, the rotation (|72[) corresponds to switching 
from Heisenberg picture to a rotating axes representa- 
tion. The generator of the unitary transformation is the 
Hcrmitian operator 

K = i{¥$> - = t(a[a 21 - a L ai 2 )/VNA- (75) 

Operators in the rotating axes representation are ob- 
tained from those in the Heisenberg picture in the fol- 
lowing way, 



-iOKa i8K 

e ze 



(76) 



The equations of motion for the rotating polariton vari- 
ables are 

^=|(Ri)=^, B i] + (iR)fi = i[4^ > (77) 



with H 



s'tHz' 



H = -AE 1 E + W(&Z- 



H = RHR + jRR 



-iG \ 
id W , (78) 
W -A 



where W = \fg 2 N A + ^ 2 . 

We immediately recognize that in the adiabatic limit 
f ? <C W, the dark-state polaritons are decoupled from 
the dynamics. Moreover, they do not involve the excited 
atomic state |3) and hence are immune to spontaneous 
emission. Therefore, if the initial state of the system con- 
sists of dark-state polaritons only and the mixing angle 
9 varies slowly enough, then the quantum state of the 
system adiabatically follows the smoothly changing dark 
eigenstates of the Hamiltonian, so the system stays in 
the same superposition of dark states as it was initially 
in. (Note that the zero eigenvalue is nondegenerate in 
each invariant subspace S n , so there is no level crossing.) 
Thus, adiabatic rotation of the mixing angle 9 from to 
7r/2 leads to a complete transfer of the photonic state to 
collective atomic excitations between levels |1) and |2). 



Indeed, for 9 = 0, the dark-state polariton mode 1 &(9) 
coincides with the signal mode cll, while it solely cor- 
responds to the collective spin excitation o^i/VA^L for 
9 = 7r/2. This is the standard way of interpreting the 
adiabatic EIT storage process [ljl [l6| . 

The bright-state polaritons are coupled to the excited- 
state mode that decays by spontaneous emission with 
rate T. If this relaxation process is fast enough (W <C T) 
or the one-photon resonance is far detuned (W <C A) , we 
can eliminate the excited mode. Then we find an effective 
decay 75 of the bright-state polaritons and a shift ljb in 
their energy, 



W 2 A 



W 2 T 



r 2 /4 + A 2 ' 



7b 



r 2 /4 + a 2 ' 



(79) 



1. Far off-resonant regime 

Now we identify two important limiting cases. In the 
first, the one-photon resonance is far detuned with re- 
spect to the Rabi frequency Vl of the control field, to the 
ensemble-enhanced vacuum Rabi frequency g\J ' Na of the 
signal field, and also to the decay rate T of the excited 
states. In this case (W <§; A and T <C A), the bright 
energy shift ojb dominates over the decay, 



r W 2 

IB = U B < OJB = • 



(80) 



Switching back to Heisenberg picture with non-rotating 
light and atomic variables, we obtain the following Hamil- 
ton operator, 



H 




(81) 



The first two terms are AC Stark shifts for the sig- 
nal mode interacting with the atoms in level |2) and 
for atoms in level |1) interacting with the intense con- 
trol field, respectively. These terms were eliminated in 
Sec. IIII Bl by considering a manifold of excited states 
and tuning the laser fields to a point where their AC 
Stark contributions exactly compensate. Finally, the 
third term in ([51]) coincides with (|53^) if we take = g\a\ 
(the factor 2 is due to the fact that there we had two 
upper levels). Finally we mention that the characteristic 
time T of the write process is determined by (|65[) and 
is of the order T - A/IF 2 . Eq. ^ then implies that 
accumulated losses are of the order 7s T ~ T/A and can 
be safely neglected. 



2. Resonant regime 

The decay of bright-state polaritons cannot be ne- 
glected near one-photon resonance. Instead, we obtain 



11 



an effective Liouvillian for the rotating density operator 
p 1 of the dark and bright polariton modes, 

^jZP 1 = -i [HeS, p] + Aiss/5'. (82) 

Here the effective Hamiltonian in the rotating axes rep- 
resentation is given by 

H cfi = oj B &$ - i9(¥$ - , (83) 
and the dissipative Lindblad superoperator is 

£& ss p=?j(2$p& -&$p-p&&). (84) 

If the adiabatic condition 9 <C Jb is satisfied, terms 
with 9 can be neglected, so the dark and bright polari- 
tons become adiabatically decoupled. This means that 
in the rotating axes representation, both the dark and 
bright polaritons can be considered as free quasi-particles 
(except that the latter have finite lifetime). If the initial 
state contains no bright polaritons at all — that is, the 
atomic ensemble is totally polarized at 9 = in the be- 
ginning of the write process (or the signal light mode 
is in the vacuum state at 9 = 5 in the beginning of 
the read-out process) — the Lindblad term does not play 
a role. As a consequence, the dynamics is entirely de- 
scribed by the Hamiltonian (|83|) that can be considered 
zero for exact one-photon resonance (A = 0). However, 
when we change from the rotating axes representation 
back to Heisenberg picture, the term neglected in the 
adiabatic approximation reappears. In the original basis 
of light and matter operators, this translates to 

cfi oft y_ i6 q J , 

and thus to 

H c s = i0(4<72i - a L CTi 2 )/ v / lV7. (85) 
Introducing the atomic and light quadrature variables as 

X A = (<T2i + ai 2 )/ V2A0i, X L = (a L + a{)/V2, 

P A = {a 21 -o 12 )li^N~ A , P L = (a L -al)/iV2, 

(86) 

we find that the effective Hamiltonian is 

H eS = 6(X A P L -P A X L ), (87) 

and we recover the ideal mapping Hamiltonian (|16p . 
Complete mapping from light to atoms or vice versa is 
achieved if the time dependence of the mixing angle is 
adjusted such that 

± l = jf ^r0(r)=0(T)-0(O). (88) 



Finally we mention that non-adiabatic losses of dark- 
state polaritons are characterized by the decay rate 

la = — « 9, (89) 

IB 

where we applied the adiabatic condition 9 <C Jb- As 
the characteristic time T of the write/read-out process is 
of the order T ~ 1/9, the accumulated losses are of the 
order Tjb ~ Jb/9 <C 1. 

IV. CONCLUSIONS 

In the present paper we have investigated quantum 
memory for light in atomic ensemble using a continuous 
variable description, i.e., position and momentum vari- 
ables for light and matter degrees of freedom. In partic- 
ular we have studied in detail two types of off-resonant 
quantum memories: one based on the quantum Faraday 
effect supplemented by measurement and feedback, and 
another involving A-type atoms and Raman scattering. 
For the first scheme we analyzed the effect of inefficiencies 
in the initial atomic state preparation and imperfections 
in the light measurement. We found that the fidelity of 
the one-pass memory scheme for storing superpositions 
of light states with very different momenta is rather low 
unless atomic state preparation can be performed with 
very high accuracy. In a two-pass scheme the necessity 
of an initial atomic spin squeezing can be avoided by 
light measurement and feedback. In a triple-pass config- 
uration, atomic state preparation as well as feedback can 
be avoided altogether. Secondly, we have investigated 
an off-resonant Raman scheme and proposed a configu- 
ration in which unwanted light shifts can be cancelled. 
Finally we have discussed near resonant quantum mem- 
ories based on electromagnetically induced transparency 
in terms of position and momentum operators of light 
and matter. We have shown that this memory can be de- 
scribed by a Hamiltonian corresponding to an ideal map 
provided that before the write and read-out processes the 
atomic ensemble or, respectively, the radiation mode are 
in the appropriate initial vacuum state. In contrast to 
the common approach involving wave equations for the 
propagating polaritons, we have developed a Hamiltonian 
formalism for the atomic and photonic quadrature vari- 
ables similar to the formalism of the first two families. 
Our results allow a straightforward comparison of the 
various continuous variable quantum memory schemes in 
the same framework. 
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